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Abstract 

Three more X particles are established in the 2014 Particle Data compared with the 2012 ones. 
There are now five established X particles named as X(3872), JT(3900), X(4260), X(4360), and 
Jf(4660). Since the first X particle Jf(3872) was discovered during the search for the remaining P 
charmonium states, it is valuable to check whether the established X particles can be explained 
by quarkonium model. In this paper, we try to calculate the mass spectra of charmonium system 
by considering spin-dependent forces deduced from one gluon exchange diagrams. The confining 
potential form is taken to be linear and the free parameters are determined by least squares method 
comparing the theoretical and the observed masses of charmonium states. 
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I. INTRODUCTION 


The analysis of charmonium system motivated important steps to understand quanti¬ 
tatively the mechanism of quark confinement in various hadrons. The first approach was 
made by introducing confining potential to estimate the masses of excited states just after 
the discovery of J/^ |lj. Because the charm quark was taken to be massive, nonrelativistic 
Schrodinger equation was applied to predict the excited states and the predictions turned 
out to be quite successful. Later the same formalism was used to calculate the detailed 
splittings of spectra in the more massive system of bottomonium. These successes lead 
to the consideration of spin splittings and a systematic derivation was possible with the 
introduction of one gluon exchange diagrams |2j. 

Although the spin splittings due to spin-spin, spin-orbit, and tensor forces have been 
calculated successfully, the form of confining potential has not been determined clearly and 
various models were tried to explain quarkonium spectra. The discrimination of different 
potential models can be carried out by comparing the spectral pattern of excited states but 
the observations of excited states with orbital angular momentum had been tedious processes 
analyzing every decaying channels between excited states. During these analyzing processes 
new states with somewhat peculiar aspects of decaying patterns have been observed and 
named as X particles [3J. 

In order to predict the energy levels of quarkonium states, we need to fix several pa¬ 
rameters such as quark mass, strong coupling constant, and potential parameter. These 
parameters are usually taken to be independent from each other, however, they are interre¬ 
lated through the exchanges of self-interacting gluons. The dynamical quark mass induced 
from the excited quarkonium states includes large gluonic contributions and the strong cou¬ 
pling constant a s varies according to the value of exchanged gluonic momentum [0], and 
the potential parameter represents just the effects of gluonic interactions. Because these 
problems are not resolved as yet from first principles, we have to fix these parameters by 
considering observed data and the most explicit data to be considered are the observed en¬ 
ergy spectra. The analysis of X particles in this viewpoint is necessary and we will carry 
out this analysis in this paper. 

In Sec|TT] we will review the observational processes of X particles, and in SecJTTTJ char¬ 
monium states are calculated with spin-dependent forces. The calculated and the observed 
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X particles are compared in SecJTV] and the final section is devoted to discussions. 


II. OBSERVED X PARTICLES 

The history of X particles is quite long. For example, in 1974 Meson Table we can find 
five states X(1430), X(1440), X(1690), X(1975), and A"(2500 — 3600). All these states 
were unestablished at that time while only 21 meson states were accepted as established 
for isoscalar or isovector states. The state X(1430) was renamed as 7/(1430) when the 
naming scheme for mesons was changed in 1986. The change of naming scheme had been 
motivated by former papers [5] classifying meson states into radially and orbitally excited 
states of quark-antiquark bound system. In doing so the spin-dependences of the observed 
meson spectra played critical roles and it turned out that the magnitudes of spin splittings 
could not be easily accounted by perturbative calculations. A typical example is the mass 
difference between 7 t( 135) and p(770) which are taken to be the spin singlet and triplet states 
of isotriplet combination of up and down quarks. This example clearly shows the problem 
of definition of quark mass generating different viewpoints between current quark mass and 
constituent one. Another example of changed name was X(1440) and this state became 
p(1450) in 1988 Particle Data. 

In 1988 Particle Data, we have another five states A"(1700), X(1850), X(1935), A"(2220), 
and X(1900 — 3600) [6]. Of these X(1850) became 0 3 (185O) and X(2220) was changed into 
fj( 2220) in 2002 Particle Data just before the discovery of X(3872) in 2003. The state 
X(1700) was replaced by X(1600) and X(1935) became X(2000) in 2002 Particle Data jTj 
but both states disappeared after the discovery of X(3872). The state X(1900 — 3600) 
existed in 1974 as A"(2500 — 3600) which had mass range overlapping with IS and IP states 
of charmonium system. The narrowness of the widths of these charmonium states prevented 
them to be checked as independent states at that time and the extension of the mass range 
from 2500 MeV to 1900 MeV could be related to the appearance of charmed particles. 
Anyway only two X particles X(1600) and X(2000) remained in 2002 Particle Data, but 
the name A" was assigned only to A"(3872) after its discovery. In 2004 Particle Data, only 
one state X(3872) was named as X particle and the meaning of the assignment X had been 
changed into the designation for non-quarkonium state (8j. 

X(3872) was discovered during the search for the 2P charmonium states and this state 
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had been established since 2006. In 2006 Particle Data [9], one component of 2P charmonium 
state \c 2 had been tabulated also and it was necessary to check whether the other components 
could be detected independently from the observed X(3872). These searches resulted in the 
report of F(3940) and H(4260). Two years later, F(3940) had been changed into A"(3940) 
and F(4260) became A"(4260) [TO] . The state E(3940) had been found in the uJ/ip invariant 
mass distribution for exclusive B decays into KuJ/ip by Belle Collaboration [IT]. But the 
state is still not established because of lack of confirmations by other experimental groups. 
On the other hand, the state A(4260) was found by BaBar Collaboration in the invariant- 
mass spectrum of 7r + 7r ~J/ip through the study of initial-state radiation events [ 12 ]. This 
state had been established in the 2008 Particle Data. 

Another state reported in 2006 Particle Data is X(1835). This state was observed by 
BES Collaboration in the invariant mass spectrum of pp pairs from radiative J/ip decays 
p~3 j. The first suggestion to account for the decay channels was a pp baryonium which 
could be a tetraquark generator. However, other possibilities such as pseudoscalar glueball 
or radial excitation of rf had been considered also but the state could not be established 
because other groups could not confirm the existence of it. Instead of confirmation of the 
same state, another state A"(1840) had been added to the 2014 Particle Data [14] . 

In 2008 Particle Data, two new state X(3945) and X(4360) had been reported and the 
state F(4260) was changed into the established A"(4260). The state X(3945) appeared in 
2010 Particle Data [TS] but disappeared after that. X(4360) stayed as unestablished state for 
several years and finally became established one in 2014 Particle Data. Similarly X(4660) 
appeared in 2010 Particle Data and became established state in 2014 Particle Data. In 
addition there appeared many X particles in 2010 Particle Data such as A"(4050), X(4140), 
X(4160), A(4250), A"(4350), and A"(4430). These particles are not established as yet but 
the importance to analyze the spectra with respect to a given model is increased. In 2012 
Particle Data jTB], three more X particles A"(3915), X(10610), and A"(10650) were added 
and the established state A"(3915) was assigned to 2 P Xco state. Finally, four more states 
A"(1840), A"(3823), A"(3900), and A"(4020) were added to 2014 Particle Data and these are 
tabulated in Table 1 according to their appearance in Particle Data. 
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Table 1. X particle states reported in Particle Data. 


(years) 







1974 

1988 

2002 

2004 2006 2008 

2010 

2012 

2014 

A(1430) 

A(1440) 


A(1600) 





A(1690) 

A(1700) 


A(1835) A(1835) 

A(1835) 

A(1835) 

A(1835) 


A(1850) 

A(1935) 





A(1840) 

A(1975) 

A(2220) 

A(2000) 





A(2500 - 3600) A(1900 - 3600) 












A(3823) 




A(3872) *A(3872) *A(3872) 

*A(3872) 

*A(3872) 

*A(3872) 

*A(3900) 






*A(3915) 





Y (3940) A(3940) 

A(3940) 

A(3940) 

A(3940) 




A(3945) 

A(3945) 


A(4020) 





A(4050) 

A(4050) 

A(4050) 





A(4140) 

A(4140) 

A(4140) 





A(4160) 

A(4160) 

A(4160) 





A(4250) 

A(4250) 

A(4250) 




F(4260) *A(4260) 

* A (4260) 

*A(4260) 

*A(4260) 





A(4350) 

A(4350) 

A(4350) 




A(4360) 

A(4360) 

A(4360) 

*A(4360) 





A (4430) 

A(4430) 

A (4430) 





A(4660) 

A(4660) 

*A(4660) 






A(10610) 

A(10610) 






A(10650) 

A(10650) 


* represents established states. 


III. CHARMONIUM STATES WITH SPIN-DEPENDENT FORCES 


The quarkonium system composed of one quark-antiquark pair with masses rn j and m 2 
can be described by the equation 




( 1 ) 
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where the Hamiltonian H can be split as 


H = H 0 + e(r) + V SD 


( 2 ) 


with 


Ho = \ m( + pf + \/mf 2 + p; 


( 3 ) 


The potential e(r) represents the spin-independent part and Vsd is the spin-dependent po¬ 
tential. By calculating relativistic propagator corrections in the Wilson loop formed by the 
quark-antiquark pair, we obtain [2J 


Vsd(t) = - 


1 f Si ■ L s 2 • L\ /I de{r ) 2 dV i 


9 

m\ 


+ 


ml 


r dr + r dr 


1 1 d \/2 

H-(si + s 2 ) • L ——— 

mim 2 r dr 


3m i m 2 

with the potentials defined by 

3 - e(r) r si • L l dVi(r) _ ig^ rT 


(3si • r s 2 ■ r - Si • s 2 )V r 3 (r) + —Si ■ s 2 H 4 (r) 


3mim 2 


( 4 ) 


Te e(r;i ? ■ xv / _ ^ ^ — t)Si • (B(x 1 ,t)E(x 1 ,t / )) • V(xi,0) , 

mf r dr J 0 Jo 

Te -dr)T s^ld^r) = - r T dt f T dt , si ( B(x ^ E ( x . v (x 2 , T) , 

mim 2 r dr mrm 2 Jo Jo 


Te~ e{r)T 


mim 2 


(Si • f s 2 • r - ~si ■ s 2 )V 3 (r) + • s 2 H 4 (r) 


.9 


f dt f dt' (Si • B(x 1; t)s 2 • B(x 2 ,t')) . 
9o 9o 


m\m 2 Jo 


The expectation values are to be evaluated along the Wilson loop as 


(6(x)) = / [dA*] Tr { P 


exp l i rj 


ig (f dz^{z)\ 6{x) 


J C 


iSym 


(6) 


( 6 ) 


where P represents the path-ordered exponential, and the non-Abelian electric and magnetic 
fields E and B are given by 


E = —VA° — A — ig [A, A°] , 
B = V x A + igA x A 


( 7 ) 


in terms of the gluonic field A 11 . The calculation of VJs can be carried out by estimating the 
gluonic correlation functions and because these correlation functions include nonperturbative 
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contributions we have to introduce some method to parametrize them. One method to 
calculate them is to use lattice QCD and there were attempts to estimate the correlations 
up to 0.6fm or so 02]' However, in order to check out whether the observed X particles 
can be assigned to quarkonium states or not it is better to use explicit form of confining 
potential and perturbative approximations to Vi’s. In this paper we will take 


r 4 1 

e ( r ) = ~2~ + b , 

a 2 3 r 

and perturbative results are 


1 dV 2 (r) _ 4 1 

r dr 3 ^r 3 ’ 


V:${r) = 77^4 , 


V 4 (r) = -a s 87iS(r) 


( 8 ) 


(9) 


The form of Vi (r) is fixed by the well-known relation [ 18J 


e(r) = V 2 (r) - V x {r) 


( 10 ) 


Then the spin-dependent potential becomes 


T/ ( \ 1 f Sl ' L , s 2 • L 

Vsd\t) = - ( —5- + 


2 V rn'i 


1 4 1\ 4 

--1- I + xQ). 

3 mi77l2 


a 2 r 3 r 3 


(si + S2) ■ L — 


4 1 . 1 

+ xcts-( 3s i ■rs 2 -r-s r S2 - + 

3 mi??T2 f 


3 1 3m^ Sl ' S 4 Mrf(r) ' (U) 


In order to solve the Eq.(l), we have to expand the square root operators in Eq.(3) by 
introducing the expansion parameter M 


M = J < p 2 > +m. 


( 12 ) 


including the momentum expectation value. M is called the effective mass and we get 


p? + ™? = f + 


1717 


+ 


PJ 


2 M 2 M 


(13) 


and the Eq.(l) can be reduced to Schrodinger-like equation. But now the singular behaviors 
in -4 potential and 5(r) have to be modihed in some way and we adopt the slight change in 
dj as ^ r+ ^, p- and the <5(r) is replaced by 

1 _7rr^ 

5(r) = -g e do" . (14) 

r o 
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For the system of charmonium states, we can take mi = m 2 and the effective masses are 
changed from state to state. However, in order to check whether the X particles can be 
assigned to charmonium states it is better to reduce the number of parameters, and so we 
will take a typical effective mass by comparing observed masses with calculated masses. The 
comparison is carried out by least squares method with variation of parameters. The values 
of r q and tq turn out to be not so much affective to the energy spectra and therefore we will 
fix these values as 0.01 and 1.0 GeV -1 . The parameter b is determined by fitting the energy 
eigenvalue of l 3 Si state to the observed mass of J/ip. Then the remaining three parameters 
are charm quark mass m, potential parameter a, and strong coupling constant a s . These 
parameters are fully varied to find the best fit to the observed charmonium spectra. For 
the fixing of parameter values, we used 10 (cc) states definitely assigned to each spin-orbit 
states of charmonium in the 2014 Particle Data. The states are listed in Table ED 


Table 2. 10 (cc) states definitely assigned to each spin-orbit states of charmonium in the 
2014 Particle Data. 


Name 

I G (J PC ) 

Mass (MeV) 

Spin-orbit states 

hc(lS') 

0+(0“ + ) 

2983.6 

1 '57 

JM1S) 

0-(l—) 

3096.916 

l 3 Si 

Xco(lP) 

0+(0 ++ ) 

3414.75 

1 3 Po 

Xcl(lP) 

0+(l++) 

3510.66 

1 3 P 1 

h c (\P) 

? ? (1+-) 

3525.38 

1 'Pi 

Xc 2 (l P) 

0+(2++) 

3556.2 

1 3 Po 

Vc(2S) 

0+(0-+) 

3639.4 

2 'S () 

<A(2 S) 

0(1 ) 

3686.109 

2 3 Si 

Xco(2P) 

0+(0 ++ ) 

3918.4 

2 3 P 0 

Xc2(2P) 

0+(2++) 

3927.2 

2 3 P 2 


The three parameters are determined by finding out the least value of 


AM 


1 

N 


E ( B -“‘ 



(15) 


where Ef al are the calculated masses of charmonium states and E° bs are the observed masses. 
N is the number of states used to calculate AM and for the fixing of parameters we have 
chosen N = 10. The variations of AM with respect to each parameter are shown in Fig.l, 
Fig.2, and Fig.3. The least value of AM is 


AM = 36.715 MeV 


(16) 
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with the parameter values 


m = 4.3 GeV, a = 2.0 GeV“\ 
With the determined parameters we can 

Table [37] 

46 

44 

¥ 42 

2 

S 

< 40 

38 

36 


Fig. 1. The variation of AM with respect to m. 



Quark mass [GeV] 


cr s = 0.28 . (17) 

calculate the masses of charmonium states as in 


> 

CD 


§ 

<1 



Fig. 2. The variation of AM with respect to a. 
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Fig. 3. The variation of AM with respect to a s . 


IV. COMPARISON OF OBSERVED STATES WITH CALCULATED RESULTS 


There are 9 more established states that can be assigned to charmoninm states if possible. 
They are listed in Table 4 with possible assignments. In order to check whether one observed 
state can be assigned to some charmonium state or not, it is valuable to calculate AM with 
inclusion of that state to the list of fixed states in Table ED Of course there should not be 
overlapping of states, and therefore -0(3770) cannot be assigned to 2 3 Si state because of 
0(2 S) and A'^OO) 3 " has to be excluded from 2 3 P 0 and 2 3 P 2 states because of x c o(2P) and 
X c2 (2P). Now we can calculate the least square difference AM for each inclusion of observed 
state and the results are shown in Table l5~~l Since the AM for the 10 fixed states is 36.715 
MeV, it seems to be appropriate to take the upper limit of 50 MeV to discriminate the 
assignments. Then 0(4160) cannot be assigned to 3 3 Si state and X(4260) is excluded from 
any assignment to charmonium states. Now except for X(3872) the remaining states have 
two possible (cc) assignments and we try to find out the best assignment by considering 
each combination of possible assignments. Firstly, the state A"(4660) can be assigned to 
5 3 Si and 4 3 Di states and we can compare these assignments by calculating AM. Because 
the state A"(3872) can be overlapped on 2 3 P X state with A"(3900) =t , we will consider these 
states later and the first comparison is carried out with only 0(3770) fixed at 1 3 Pi state. 


The results are shown in the first two rows of Table 6. It turns out that A"(4660) is better 
suited to 5 3 Si than 4 3 Pi, so we can fix X(4660) as 5 3 Si (cc) state. The next comparison 
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Table 3. The caculated masses of cc states with m = 4.3 GeV, a = 2.0 GeV ', a s = 0.28 . 


n 

states 

Energy 

(MeV) 

n 

states 

Energy 

(MeV) 

n 

states 

Energy 

(MeV) 

n 

states 

Energy 

(MeV) 

i 

'So 

3069.289 

2 

3 Po 

3861.784 

i 

3 d x 

3758.536 

i 

3 Fs 

4013.832 

2 

'So 

3658.608 

3 

3 Po 

4211.973 

2 

3 D\ 

4119.056 

2 

3 Fs 

4341.675 

3 

'So 

4063.475 

4 

3 Po 

4522.449 

3 

3 D\ 

4436.514 

3 

3 Fs 

4638.590 

4 

'So 

4402.871 

5 

3 Po 

4806.507 

4 

3 D x 

4725.683 

1 

3 F 4 

4054.074 

5 

'So 

4705.338 

1 

3 Pi 

3502.678 

1 

3 d 2 

3787.032 

2 

3 f 4 

4382.869 

6 

'So 

4983.034 

2 

3 P1 

3921.409 

2 

3 d 2 

4148.184 

3 

3 f± 

4680.448 

1 

3 Si 

3096.916 

3 

3 Pi 

4269.803 

3 

3 d 2 

4466.071 

1 

'G 4 

4223.334 

2 

3 S! 

3670.321 

4 

3 Pi 

4578.827 

4 

3 d 2 

4755.562 

2 

'g 4 

4528.210 

3 

3 S! 

4071.601 

5 

3 Pi 

4861.634 

1 

3 d 3 

3823.186 

1 

3 g 3 

4175.334 

4 

3 S! 

4409.295 

1 

3 p 2 

3549.774 

2 

3 d 3 

4185.077 

2 

3 g 3 

4479.266 

5 

3 Si 

4710.737 

2 

3 p 2 

3966.525 

3 

3 d 3 

4503.443 

1 

3 G4 

4214.091 

6 

3 Si 

4987.737 

3 

3 P 2 

4314.022 

4 

3 d 3 

4793.273 

2 

3 g 4 

4518.777 

1 

1 P1 

3526.129 

4 

3 P 2 

4622.503 

1 

1 P's 

4023.559 

1 

3 g 5 

4261.390 

2 

1 PI 

3944.129 

5 

3 P 2 

4904.928 

2 

1 P's 

4351.657 

2 

3 g 5 

4566.997 

3 

1 PI 

4292.203 

1 

1 d 2 

3798.540 

3 

'P'S 

4648.749 

1 

'G 5 

4406.616 

4 

1 PI 

4601.036 

2 

'£>2 

4159.998 

1 

3 f 2 

3981.946 

1 

3 g 4 

4350.730 

5 

1 P1 

4883.709 

3 

'£>2 

4478.081 

2 

3 f 2 

4309.036 

1 

3 g 5 

4397.523 

1 

3 Po 

3440.481 

4 

'£>2 

4767.712 

3 

3 f 2 

4605.420 

1 

3 g 6 

4452.983 


Table 4. 9 states that can be assigned to cc states. 


Name 

I G (J PC ) 

Mass(MeV) 

Possible (cc) states 

0(3770) 

o-(i-) 

3773.15 

2 3 S'i , 1 3 D\ 

A(3872) 

0+(l++) 

3871.69 

2 3 P 1 

A^OO)* 

?(1+) 

3888.7 

2 1 P 1 , 2 3 P 0 , 2 3 Pi, 2 3 P 2 

0(4040) 

o (i ) 

4039 ± 1 

2 3 D\ , 3 3 Si 

0(4160) 

o-(i-) 

4191 ± 5 

2 3 D\ , 3 3 Si 

A (4260) 

? ? (1--) 

4251 ± 9 

2 3 Di , 3 3 Si, 3 3 P»i, 4 3 5i 

A (4360) 

? ? (1--) 

4361 ± 13 

3 3 D\ , 4 3 S\ 

0(4415) 

o-(i-) 

4421 ± 4 

3 3 P>i , 4 3 Si 

A(4660) 

? ? (1-) 

4664 ± 12 

4 3 D\ , 5 3 Si 


is the assignments of 3 3 Si and 2 3 D 3 states to "0(4040) and 0(4160). These two cases 
are compared in the third and the forth rows of Table 6. The case of 3 3 S'i to 0(4160) 
and 2 3 D 4 to 0(4040) is excluded because AM exceeds 50 MeV. Therefore we can now fix 
0(4040) to 3 3 S'i state and 0(4160) to 2 3 D\ state. Another comparison is between the 
assignments of 4 3 S 1 and 3 3 D 1 states to AT(4360) and 0(4415). These cases are classified 


10 
























Table 5. Calculation of AM of 9 states that be assigned to cc states and assignment 
possibility 


Name 

Assigned state 

AM(MeV) 

Possibility 

V>(3770) 

1 3 Di 

35.283 

o 

A (3872) 

2 3 P 1 

38.082 

o 

A(3900) ± 

2 l P 1 

38.792 

o 

2 3 Pi 

36.370 

o 

-0(4040) 

3 3 Si 

36.361 

o 

2 3 Di 

42.522 

o 

0(4160) 

3 3 Si 

50.214 

X 

2 3 Di 

41.183 

o 

A (4260) 

3 3 Si 

64.431 

X 

4 3 Si 

59.190 

X 

2 3 Di 

52.992 

X 

3 3 Di 

65.986 

X 

A(4360) 

4 3 SI 

37.915 

o 

3 3 Di 

41.760 

o 

A(4415) 

4 3 Si 

35.184 

o 

3 3 Di 

35.318 

o 

A(4660) 

5 3 Si 

37.737 

o 

4 3 Di 

39.640 

o 


Table 6. Combination of possible assignments and the calculation of AM 


List 

Charmonium states 

AM 

3 3 Si 

4 3 51 

5 3 5i 

2 1 P 1 

2 3 Pi 

1 3 D 1 

2 3 D 1 

3 3 L>i 

4 3 D 1 

(MeV) 

1 



A (4660) 



0(3770) 




36.376 

2 






0(3770) 



A (4660) 

38.187 

3 

0(4040) 


A (4660) 



0(3770) 

0(4160) 



39.746 

4 

0(4160) 


A (4660) 



0(3770) 

0(4040) 



51.090 

5 

0(4040) 

0(4415) 

A (4660) 



0(3770) 

0(4160) 

A(4360) 


41.800 

6 

0(4040) 

0(4415) 

A (4660) 


A (3872) 

0(3770) 

0(4160) 

A(4360) 


42.307 

7 

0(4040) 

0(4415) 

A (4660) 

A(3900) ± 

A (3872) 

0(3770) 

0(4160) 

A(4360) 


43.141 

8 

0(4040) 

0(4415) 

A (4660) 


A(3900) ± 

0(3770) 

0(4160) 

A(4360) 


41.321 

9 

0(4040) 

A(4360) 

A (4660) 



0(3770) 

0(4160) 

0(4415) 


39.283 

10 

0(4040) 

A (4360) 

A (4660) 


A(3872) 

0(3770) 

0(4160) 

0(4415) 


39.972 

11 

0(4040) 

A(4360) 

A (4660) 

A(3900) ± 

A (3872) 

0(3770) 

0(4160) 

0(4415) 


40.984 

12 

0(4040) 

A (4360) 

A (4660) 


A(3900) ± 

0(3770) 

0(4160) 

0(4415) 


38.927 


according to the remaining assignments of 2 1 P i and 2 3 Pi states to A^OO^ and A(3872). 


They are shown in the lower part of Table 6. It turns out that A"(4360) is better suited 
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to 4 3 S i state and -0(4415) is naturally assigned to 3 3 D\ state. With these assignments 
the best fit to the observed states has AM = 38.927 MeV not so different from the value 
AM = 36.715 MeV with only 10 fixed states. The best fit includes A"(3900) =b as 2 3 Pi state 
and therefore A(3872) cannot be included in the spectra of charmonium states. 

We analyzed only five established X particles. Of these five X particles, two X particle 
A"(3872) and A(4260) turn out to be non-(cc) states. However, there exist 18 X particles in 
2014 Particle Data and it seems impossible to explain them with quarkonium model. Other 
possibilities such as hybrid states or tetraquark states have to be considered to account for 
these observed states. In order to check whether the other X particles can be assigned 
to charmonium states or not, we need to repeat the above processes. But for the other 
unestablished X particles it may be sufficient to compare with the calculated spectra like 
in Table HD Because we have obtained the best fit to the observed established states, it 
is necessary to recalculate the whole spectra with new parameters determined from the 17 
states assigned to charmonium states. The new parameters are 

m = 5.0 GeV, a = 2.0 GeV~\ a s = 0.28, (18) 

and the calculated results are shown in Table EH 

V. DISCUSSIONS 

In this paper, we have calculated the energy spectra of charmonium system by determin¬ 
ing the parameters with least squares method. The determined quark mass turns out to be 
very large compared wirh the current quark mass deduced from the electroweak interactions. 
The difference results from gluonic interactions that increase with higher excited energies. 
The effect of large dynamical mass appears as reduction of spin splittings as they are propor¬ 
tional to the inverse square of quark mass. Thus the largest difference between the calculated 
and the observed masses of rj c (lS) can be ascribed to the large quark mass determined to fit 
all the observed data. Moreover, for S-wave states, the probability of quark pair annihilation 
is not negligible and similar considerations lead to the introduction of vacuum condensate 
idea to predict the splittings between the triplet and the singlet states. The second state 
with large difference between the calculated and the observed masses is 0(4160). This state 
is assigned to 2 3 D 1 and it seems necessary to check the change of measured mass from the 
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Table 7. The calculated results with m = 5.0 GeV, a = 2.0 GeV 1 , a s = 0.28. 


States 

Theory 

(MeV) 

cc 

States 

Theory 

(MeV) 

CC 

Name 

Mass (MeV) 

I G (J PC ) 

Name 

Mass(MeV) 

I G {J PG ) 

1 % 

3072.704 

^c(lS) 

2983.6 

O ' (0 " G 

2 1 D 2 

4144.799 




1 3 S 1 

3096.916 

JMis) 

3096.916 

o-(i — ) 

2 3 Pi 

4112.396 

^(4160) 

4191±5 

0"(1“) 

2 % 

3655.824 

r? c (2 S) 

3639.4 

0+(0-+) 

2 3 D 2 

4135.587 




2 3 Si 

3665.211 

V>(2 S) 

3686.109 

o-(i —) 

2 3 D 3 

4164.523 




3% 

4047.043 




3 1 D 2 

4449.003 




3 3 Si 

4053.409 

^(4040) 

4039±1 

o-(i-) 

3 3 P>i 

4416.202 

V>(4415) 

4421±4 

o-(i-) 

4 % 

4372.978 




3 3 D 2 

4439.656 




4 3 Si 

4377.961 

X(4360) 

4361T13 

? ? (1--) 

3 3 D 3 

4468.898 




5 % 

4662.635 




4 1 D 2 

4725.639 




5 3 S ± 

4666.803 

X(4660) 

4664±12 

? ? (1--) 

4 3 D\ 

4692.543 




6 % 

4928.151 




4 3 D 2 

4716.197 




6 3 Si 

4931.773 




4 3 D 3 

4745.657 




1 l Pi 

3532.609 

h c (\P) 

3525.38 

? ? (1 + -) 

1 3 p 3 

4015.570 




1 3 Po 

3459.000 

Xco(l P) 

3414.75 

0+(0++) 

1 3 p 2 

3983.757 




1 3 P 1 

3512.699 

Xcl(l P) 

3510.66 

0+(0++) 

1 3 p 3 

4008.207 




1 3 p 2 

3552.889 

Xc2(lP) 

3556.2 

0+(2++) 

1 3 P 4 

4038.824 




2 : Pi 

3935.717 




2 3 P 3 

4329.260 




2 3 P 0 

3866.116 

Xco(2P) 

3918.4 

0+(0++) 

2 3 F 2 

4296.696 




2 3 Pi 

3916.802 

x^o)* 

3888.7 


2 3 P 3 

4321.706 




2 3 P 2 

3954.587 

Xc 2 (2P) 

3927.2 

0+(2++) 

2 3 P 4 

4353.035 




3 : Pi 

4269.587 




3 3 P 3 

4612.958 




3 3 P 0 

4202.176 




3 3 F 2 

4579.867 




3 3 Pi 

4251.117 




3 3 F 3 

4605.269 




3 3 P 2 

4287.814 




3 3 P 4 

4637.096 




4 ^ 

4565.092 




1 3 G 4 

4207.537 




4 3 P 0 

4499.221 




1 3 g 3 

4171.471 




4 3 Pi 

4546.888 




1 3 G 4 

4200.635 




4 3 P 2 

4582.929 




1 3 g 5 

4236.089 




5 : Pi 

4835.191 




2 3 G 4 

4498.623 




5 3 P 0 

4770.541 




2 3 G 3 

4461.832 




5 3 Pi 

4817.173 




2 3 G 4 

4491.573 




5 3 P 2 

4852.756 




2 3 G 5 

4527.735 




1 1 d 2 

3798.327 




1 

4383.196 




1 3 £>i 

3766.507 

^(3770) 

3773.15 

o-(i —) 

1 3 if 4 

4341.549 




1 3 p 2 

3789.320 




1 3 #5 

4376.447 




1 3 d 3 

3817.792 




1 3 ii 6 

4417.722 
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value around 4160 MeV. 

The main conclusion of this paper is that the two established X particle X(3872) and 
A"(4260) cannot be assigned to charmonium states. As is well known from the time of 
discovery, A^(3872) is considered to be a state of tetraquark. However, the calculation of 
energy states for 4-quark degree of freedom is not an easy one and the observed mass of 
X(3872) is not quite different from the mass of 2 3 Pi charmonium state. In 2014 Parti¬ 
cle Data, A"(3900) ± state is established and stays nearer to 2 3 Pi charmonium state than 
A"(3872), but it is possible for A"(3900) ± to be assigned to 2 1 P 1 state and then A"(3872) 
can be assigned to 2 3 Pi charmonium state. Further study is needed by considering decay 
processes and state mixings to clarify the status of A"(3872). In contrast, AT(4260) is clearly 
excluded from the assignments to charmonium states and therefore it is plausible to check 
the calculations of tetraquark states to accommodate X(4260). 

Another state to note is x c o(2P). This state was named as A"(3915) in 2012 Particle 
Data and established immediately after its discovery. However, the difference between the 
observed mass and the calculated mass in Table 3 amounts to 56 MeV. It may be possible 
to assign X(3915) to 2 3 Po charmonium state but the significant difference between the 
calculated and the observed masses could be taken to indicate other possibilities. Other X 
particles not discussed in this paper are still unestablished and we need to confirm these 
states with more experiments. Peculiar examples are A’(1835) and A(1840). These states 
cannot be assigned to any charmonium state and it is open to question whether they can be 
accounted as tetraquark state or as even more complex combination of mesons and baryons. 
Two heavy X particles X(10610) and X(10650) are thought to be some states containing 
bottom quark and we need to analyze bottomonium system. 

In summary, now there exist so many X particles that they cannot be accommodated 
as quarkonium sates. Of course some X particles can be explained by quarkonium states 
but we need to introduce other possibilities such as tetraquark states, hybrids including 
gluonic degrees of freedom, and so on. In order to improve theoretical calculations, we need 
to rederive the form of gluon propagator in space-time coordinates to get rid of the diver¬ 
gence behavior at short distances. More systematic derivations of potential forms are also 
long-standing problems in strong interaction phenomena. Further establishments of more 
X particles may generate active researches on such subjects resulting in the quantitative 
understanding of strongly interacting bound systems. 
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